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Abstract. It is shown that many properties of Lya forest absorbers can be de¬ 
rived using simple physical arguments. Analytical expressions are derived for the 
density and the size of an absorber as a function of its neutral hydrogen column den¬ 
sity, which agree well with both observations and hydrodynamical simulations. An 
expression is presented to compute flioM from the observed column density distribu¬ 
tion, independent of the overall shape of the absorbers. Application to the observed 
column density distribution shows that at high redshift most of the baryons are in 
the forest and suggests a simple interpretation for its shape and evolution. 


1 Introduction 

In the last decade semi-analytic models (e.g., [||) and hydrodynamical simula¬ 
tions (see [|j for a recent review) have been used to show that cold dark matter 
models with a nearly scale-invariant spectrum of initial, adiabatic, Gaussian 
fluctuations are very successful in reproducing the large body of high-quality 
observations of the Lycr forest. The physical picture that has emerged is that 
the forest arises in a network of sheets, filaments, and halos, which give rise 
to absorption lines of progressively higher column densities. The low column 
density absorption lines {Nhi cm“^ at z ~ 3), arise in a smoothly 

fluctuating intergalactic medium (IGM) of moderate overdensity (p < 10 (p)), 
and contain most of the baryons in the universe. On large scales the gas traces 
the dark matter and observations of the Lya forest can be used to study the 
large-scale distribution of matter, while on small scales hydrodynamical ef¬ 
fects are important and the detailed line profiles can be used to reconstruct 
the thermal history of the IGM. In this contribution I will show that the mod¬ 
ern picture of the forest can be derived directly from the observations using 
straightforward physical arguments, without making any assumptions about 
the presence of dark matter, the mechanism for structure formation or the 
precise cosmological model. This work is described in more detail in a recent 
publication [||. 


2 Physics 

Gonsider an (optically thin) self-gravitating gas cloud of arbitrary shape that is 
intersected by our line of sight to a background QSO. Let nn be the character¬ 
istic total hydrogen number density of the absorber, i.e., the density weighted 


by the neutral hydrogen column density and let L be the characteristic size of 
the absorber, i.e., the size over which the density is of order uh- 

Regardless of whether the cloud as a whole is in dynamical equilibrium, 
it will in general not be far from local hydrostatic equilibrium, i.e., L ~ Lj 
where Lj is the local Jeans length. If L <C ij, then the cloud will expand or 
evaporate and equilibrium will be restored on a sound-crossing time scale. If 
L ^ Lj, then the cloud is Jeans unstable and will fragment or shock and equi¬ 
librium will be restored on a dynamical time scale. Note that this argument 
breaks down if the relevant time scale is larger than the Hubble time, because 
in that case there simply has not been sufficient time to establish equilibrium. 
Since tdyn = and tn = H~^ (G this is true for underdense 

absorbers. Such absorbers have sizes greater than the sound horizon and can 
be regarded as a “fluctuating Gunn-Peterson effect”. Systematic departures 
from local hydrostatic equilibrium (L Lj) also occur in clouds confined 
by external pressure (but the large cloud sizes derived from observations of 
quasar pairs rule out the possibility that a significant fraction of the forest 
arises in pressure confined clouds 0 ), and in rotationally supported clouds 
whose spin axis is nearly perpendicular to our sight line. 


3 Results 


Combining the condition of local hydrostatic equilibrium, L ~ Lj, with the 
expression for the neutral fraction in an optically thin plasma, nni/nH ~ 
0A6nHT^^ '^^T ^2 (where Td = T/IO^ K is the temperature and ri 2 = r/10“^^ s“^ 
is the hydrogen photoionization rate), yields expressions for the density and 
size of an absorber as a function of its neutral hydrogen column density 
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where fg w Qh/flm is the fraction of the mass in gas (excluding stars and 
molecules). From the line widths of the absorption lines it is known that the 
temperature is roughly consistent with that expected from photoionization 
heating, T 4 ~ 1. Studies of the proximity effect estimate ri 2 ~ 1 at z ~ 3 and 
ri 2 ~ 10 “^ at z ~ 0 (e.g., § 1 ). 

The density - column density relation agrees with published results from 
hydrodynamical simulations to within a factor of a few, about as well as dif¬ 
ferent simulations agree with each other. The sizes of the absorbers have not 
been investigated in detail with simulations, but equation does agree with 
the sizes derived from observations of multiple sight lines (e.g., Q). 

Using equation to compute the neutral fraction as a function of column 
density, we can compute UigM) the contribution of the Lya forest to the cosmic 







baryon density, directly from the observed column density distribution: 
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Figure la shows the observed column density distribution at z ^ 3 and the 
single power law fit of Hu et al. (1995) (dashed line). Integration of this column 
density distribution according to equation^, using (r 4 , ri 2 ,/g,/i, Ha) = 
(2.0,1.0,0.16,0.65,0.3,0.7), yields « 0.015 for loglV^j = 13.0-17.2, 

confirming that at high z most of the baryons are in the diffuse IGM. 

Figure lb shows the derived mass distribution, i.e., dHiGM/rflog(l + i5). 
The shape of the mass distribution clearly reflects the deviations from a single 
power law column density distribution and agrees reasonably well with the 
actual mass distribution in a hydrodynamical simulation (solid line), kindly 
provided by T. Theuns. The agreement becomes even more impressive if 
we compare with the distribution of gas that has a temperature logT < 4.8 
(dashed line), low enough for collisional ionization to be ineffective. This good 
agreement suggests that the observed structure in the IGM formed through 
gravitational instability in an expanding universe, as is the case in the simu¬ 
lation. 

The shape of the mass distribution derived from the observations also sug¬ 
gests a simple physical interpretation for the deviations from the single power 
law in the observed column density distribution, whose origin was not yet un¬ 
derstood. The steepening at Nhi ~ cm“^ reflects the fall-off in the 

density distribution due to the onset of rapid, non-linear collapse at i5 ~ 10. 
Finally, the flattening at Nhi ~ 10^® cm“^ can be attributed to the flatten¬ 
ing of the density distribution at 5 ^ 10^ due to the virialization of collapsed 
matter. 

The evolution of the shape of the column density distribution can now 
also be explained. From z ~ 3 to z ^ 0, the slope of the column density 
distribution is found to steepen for Nhi ~ lO^^-lO^"^ cm“^ and to flatten for 
Nhi ~ lO^^-lO^® cm“^ (e.g., §)• If gravitational instability is responsible for 
the shape of the gas density distribution, then this distribution will be similar 
at z ~ 0 and z ~ 3, provided it is expressed as a function of the density relative 
to the cosmic mean. Hence, the change in the slope of the column density 
distribution must mainly be due to the evolution of the function Nhi{S), 
which is determined by the expansion of the universe and the evolution of 
the ionizing background. Equation 0 shows that Nhi oc (1 -I- z)®/^/r. The 
neutral hydrogen column density corresponding to a fixed density contrast is 
thus about a factor 5 x 10^r(z = 0)/F(z = 3) ~ 5-50 lower at z 0 than at 
z ~ 3, and this can account for the evolution in the observed distribution. 

This work demonstrates that the main aspects of the physical picture of 
the Lya forest, including the facts that the absorbers are extended, of low 
overdensity, and contain a large fraction of the baryons, can all be derived di¬ 
rectly from the observations using straightforward analytical arguments, with- 
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Figure 1: (a) Observed column density distribution function at z ~ 3. (b) 
The mass distribution as a function of overdensity, derived form the observed 
column density distribution using equation 4 and (T 4 , ri 2 ,/g, h, Slmj = 
(2.0,1.0,0.16,0.65,0.3,0.7). The curves, which have been normalized to the 
maximum of the Hu et al. data points, indicate the mass distribution in a 
hydrodynamical simulation with the same cosmological parameters as were 
used to calculate the data points. 


out making assumptions about the presence of dark matter, the mechanism for 
structure formation, or the precise cosmology. Although this analysis is essen¬ 
tially model-independent, the shape of the derived mass distribution strongly 
suggests that gravitational instability in an expanding universe is responsible 
for the distribution of matter in the IGM, thus lending credence to the use of 
ab-initio hydrodynamical simulations to investigate the detailed properties of 
the forest and their dependence on the parameters of the models. 
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